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Abstract 

Using the modified Prelle- Singer approach, we point out that exphcit time inde- 
pendent first integrals can be identified for the damped linear harmonic oscillator 
in different parameter regimes. Using these constants of motion, an appropriate 
Lagrangian and Hamiltonian formalism is developed and the resultant canonical 
equations are shown to lead to the standard dynamical description. Suitable canon- 
ical transformations to standard Hamiltonian forms are also obtained. It is also 
shown that a possible quantum mechanical description can be developed either in 
the coordinate or momentum representations using the Hamiltonian forms. 
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1 Introduction 



In a recent paper [1,2], we have studied the dynamical aspects of the Lienard 
type nonlinear oscillator equation 



X + kxx H x^ + Ax = 0, (1) 

9 

where over dot denotes differentiation with respect to t and k and A are arbi- 
trary parameters, and shown that the frequency of oscillations of a nonlinear 
oscillator need not depend upon the amplitude [1]. Our studies have also re- 
vealed that the system ([T]), eventhough deceptively appears to be a dissipative 
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one, admits conservative Hamiltonian description for all values of the parame- 
ters k and A. In fact, the conventional criterion that the divergence of the flow 
function vanishes for a Hamiltonian system fails here in the original variables. 
These facts motivate one to investigate further the possible existence of time 
independent Hamiltonian functions in actual dissipative systems. 

For our present study, we consider the ubiquitous damped linear harmonic 
oscillator 



X = —{ax + Xx) = x), (2) 

where a and A are the damping parameter and strength of the forcing, re- 
spectively. When damping is absent, a = 0, the frequency of oscillation 
ui = ^/X > 0. For A < 0, one has unbounded motion. System ([2]) is known 
to possess a time- dependent Lagrangian [3-5,7,15], L = e"*(|x^ — ^x"^), and 
a Hamiltonian H = |p^e~"* + |a:^e"*. For the past several decades there has 
been a number of attempts to quantize the damped linear harmonic oscillator 
(see for example, Refs. [6-20]) from different points of view, but it appears 
that the problem still eludes a completely satisfactory resolution. One of the 
major conceptual difficulties seem to be the lack of an appropriate time in- 
dependent Hamiltonian formalism. In this paper, we report such a formalism 
by making use of the recently developed modified Prelle-Singer approach to 
identify integrals of motion and integrability of dynamical systems [21]. 

System ([2]) is obviously a dissipative one for a > (and growing one for 
a < 0), since ^(|i^ + |Ax^) = ^ = -ax'^ < 0, where E = {^x"^ + |Aa:^) is 
the total energy of the undamped oscillator. From a dynamical point of view 
Eq. ([2]) can be written as a system of coupled first order ordinary differential 
equations (ODEs) 



x = y (3) 
y = -Xx - ay (4) 

so that the divergence of the fiow function 



A=* + *=-a<0. (5) 

OX ay 

which is the standard criterion for dissipative systems. Yet, we show here that 
system ([2]) possesses a time independent integral of motion for all values of 
the parameters A and a. We also provide the explicit form of the Lagrangian 
and Hamiltonian associated with Eq. ([2]) in different parametric regimes and 
derive all the expected solutions from the resulting canonical equations of 
motion. Our results also show that the dissipative dynamical systems can 
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also have conservative Hamiltonian description, which can aid in appropriate 
quantization procedure as will be indicated. 

The plan of the paper is as follows. In Section [21 we point out how the time 
independent integrals can be obtained for the damped linear harmonic os- 
cillator ([2]) . Then in Section [31 we identify an appropriate time- independent 
Lagrangian and so a time-independent Hamiltonian in the appropriate pa- 
rameter ranges. In Section [H we show that through the resulting canonical 
equations of motion, we obtain the correct explicit solutions. In Section [5l we 
succeed to identify suitable canonical transformations to identify 'standard' 
Hamiltonian description. Finally, in Section [HI we discuss briefly the possible 
quantum formulation of the problem. Section [7] summarizes our results. 



2 Extended Prelle-Singer procedure 

To begin with we assume that the ODE ([2|) admits a first integral I{t, x, x) = 
C, with C constant on the solutions, so that the total differential becomes 

dl = Itdt + Ixdx + I±dx = 0, (6) 

where each subscript denotes partial differentiation with respect to that vari- 
able. Rewriting Eq. ([2D in the form (pdt — dx = and adding a null term 
S{t, X, x)x dt — S{t, X, x)dx to the latter, we obtain that on the solutions the 
1-form 

(0 + Sx)dt — Sdx — = 0, (j) = —{ax + Xx). (7) 

Hence, on the solutions, the 1-forms ([2D and ^ must be proportional. Multi- 
plying ^ by the factor R{t, x, x) which acts as the integrating factors for ([71), 
we have on the solutions that 

dl = R{(p + Sx)dt - RSdx - Rdx = 0. (8) 
Comparing Eq. ([6D with ([HD we have, on the solutions, the relations 

It = R{<j) + xS), h = -RS, h = -R. (9) 

Then the compatibility conditions, Itx = Ixt, hx = ht, Ixx = hx, between the 
Eqs. ([ID, provide us [22] 

St + xSx + (i)Si: = -(i)x + 4>xS + s^, (10) 
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Rt + xR, + 4>R± = -{4>x + S)R, (11) 
Rx — SRx — RSx = 0. (12) 



Solving Eqs. (fT0|) - (fT2|) one can obtain expressions for S and R. It may be noted 
that any set of special solutions {S, R) is sufficient for our purpose. Once these 
forms are determined the integral of motion I{t, x, x) can be deduced from the 
relation 



/ = n - r2 - 
where 



dx, (13) 



ri = y R{(j) + xS)dt, r2 = J (RS + ^n)dx. 

Equation (fT3|) can be derived straightforwardly by integrating the Eq. (Q. 

As our motivation is to explore time independent integral of motion for the 
Eq. ([2]) let us choose It = 0. In this case one can easily fix the null form S 
from the first equation in as 

^ ^ -0 ^ {ax + Xx) 

X X 

Substituting this form of 5* into ffTTl) we get 



\x 

xRx — (ax + \x)Rx = R- (15) 

X 



Equation (fT5|) is a first order linear partial differential equation with variable 
coefficients. As we noted earlier any particular solution is sufficient to construct 
an integral of motion (along with the function 5*). To seek a particular solution 
for R one can make a suitable ansatz instead of looking for the general solution. 
We assume R to be of the form, 

{A{x)+B{x)x + C{x)x^Y' ^^^^ 



where A, B and C are functions of their arguments, and r is a constant which 
are all to be determined. We demand the above form of ansatz [TB^] . which is 
very important to derive the Hamiltonian structure associated with the given 
equation, due to the following reason. To deduce the first integral / we assume 
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a rational form for /, that is, I = i^^, where / and g are arbitrary functions 
of X and x and are independent of t, from which we get = ^""^"^J^"^ and I± = 
f-9-j9- . From (191) one can see that R = h = lM^f9± 5" = Is 
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RS 



s. = Sfj9^ and 

X J ±9 jQx 

Ir, SO that the denominator of the function S should be the numerator 



of the function R. Since the denominater of S" is i: (vide Eq. (HM ) we fixed 
the numerator of R as x. To seek a suitable function in the denominator 



However, it is difficult 

n(x,x) ' 



initially one can consider an arbitrary form R 
to proceed with this choice of h. So let us assume that h{x, x) is a function 
which is polynomial in x. To begin with let us consider the case where h is 
quadratic in x, that is, h = A{x) + B{x)x + C{x)x'^. Since R is in rational 
form while taking differentiation or integration the form of the denominator 
remains same but the power of the denominator decreases or increases by a 
unit order from that of the initial one. So instead of considering h to be of 
the form h = A{x) + B{x)x + C{x)x'^, one may consider a more general form 
h = {A{x) + B{x)x + C{x)x'^Y, where r is a constant to be determined. The 
parameter r plays an important role, as we see below. 



Substituting f|T6|) into f|T5l) and solving the resultant equations, we arrive at 
the relation 



x{A^ + B^x + C^x^) - {ax + \x){B + 2Cx) 



-a{A + Bx + Cx^). 



(17) 



Solving Eq. (fT71) . we can fix the forms of A, B, C and r and substituting 
them into Eq. (fT6l) we can get the integrating factor R. Doing so, we find 



R = < 



X 



(Ax^ + axx + i^) 

X 



[x + 



(r-l) 



X, 



axY 



<4A 
> 4A 

a = 



where 



a 
2A 



a 



± Va^ - 4A 



(19) 



One can easily check the functions S and R given in (IT^ and (fT8|) . respectively, 
satisfy (fT2!) also. Finally, substituting R and S into the form (fT3|) for the 
integral we get 
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i log(x^ + axx + Ax^) + ^tan ^ 



ax+2Xx 



2uix 



o? < 4A (A, a 7^ 0) (underdamped) 
{^(x + fx)(x+(^ax)(i-^), 

> 4A (A, a 7^ 0) (overdamped) 
T— f — r - log(x + ^ax), 

{x+-^ax) 2 

q;^ = 4A (A, a 7^ 0) (critically damped) 



X + ax, 



x^ + Ax^, 



A = (a 7^ 0) (pure damping) 



a = (A 7^ 0) (no damping) 



(20) 



where u> = \ J{A\ — a^). One can easily check (using Eq. ([2])) that ^ = for 



each of the cases in Eq. ( 1201) . In the above, from a physical point of view we 
consider only the principle branches to avoid multivaluedness. 

Eq. fl2U]) demonstrates clearly that the damped linear harmonic oscillator ([2]) 
admits a time independent integral of motion for all values of the parameters 
a and A. 



Now, due to the fact that the integrals of motion fl20|) are time independent 
ones we can look for a Hamiltonian description for the respective equations of 
motion. In fact, in the following we do obtain the explicit time independent 
Hamiltonians for all the above cases. 



3 Hamiltonian description of the damped linear harmonic oscilla- 
tor 

Assuming the existence of a Hamiltonian 

I{x,x) = H{x,p) = px — L{x,x), (21) 

where L{x, x) is the Lagrangian and p is the canonically conjugate momentum, 
we have 

dl dH dp . ^ dL dp . 

dx dx dx dx dx 
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From (l22l) we identify 



p = [ —dx + fix) 
J X 



(23) 



where f{x) is an arbitrary function of x. Equation fl23|) has also been derived 
recently by a different methodology [23]. Hereafter, without loss of general- 
ity we take f{x) = 0. Substituting the known expression of / into (1251) and 
integrating it we can obtain the expression for the canonical momentum p. 
Substituting back the latter into the Eq. (12T|) and simplifying the resultant 
equation we arrive at the following Lagrangian 



L 



^ ( tan 



2x+ax 
2ll)x 



—) + atan~^ 



ax+2Xx 



2lux 



log(i;^ + axx + Ax^), 



r 



(2-r) 

log(i; + ^ax), 

X log(i;) — X — ax, 

x^ — Xx"^, 



and Hamiltonian 



a2 < 4A 

„2 



a2>4A (r^O,l,^^, 



a 



4A (r 



A = (r = 1) 
a = (r = 0) 



H 



^log[x'^sec^{u!xp)] — < 4A 

('■-2) 



(r-2) yt^) 

log(p) - \axp, 
+ ax, 



,2^ '2 

where p is the canonically conjugate momentum 



> 4A (r ^ 0,1,2) 
a2 = 4A (r = 2) 
A = (r = 1) 
a = (r = 0) 



(25) 



tan 



2x-\-ax 
2ujx 



< 4A 



p=}{x+ ^axY'-^\ a' > 4A (r ^ 0, 1), 
log(i), A = (r = 1) 

p = X, a = (r = 0) 

respectively, for the damped harmonic oscillator (P 



(26) 
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One can easily check that the canonical equations of motion for the above 
Hamiltonians are nothing but the equation of motion ([2]) of the damped lin- 
ear harmonic oscillator in the appropriate parametric regimes (see below for 
details). 



4 Canonical Equations and Solutions 



We now write down the canonical equations of motion obtained from the 
above Hamiltonian expressions for each of the cases separately and obtain the 
corresponding solutions, as it should be, to show the correctness of our results. 



4.1 Case: 1 < AX: 



Now from the form of the Hamiltonian given in Eq. (1251) . we first derive the 
Hamilton's equation for the case < 4A: 



X = -x{2lj tan{LJxp) — a), (27) 

1 1 

p = -p(a — 2u taniuxp)) , (28) 

2 X 



where u = |y(4A — a^). Rewriting Eq. (1271) for p and substituting it into 
Eq. (l28l) . one can immediately check that the equation of motion ([2]) indeed 
follows straightforwordly from the canonical Eqs. ( |271) - (!28l) . Now multiplying 
Eq. fl28l) by x and fl27|) by p and adding and integrating the resultant equation 
we get 



P=^, (29) 

X 

where 6 is an integration constant. Substituting the expression for p into 
Eq. (127|) we obtain the following equation for x, namely. 



X = -x{2ujta.n{uj{6 — t)) — a), (30) 

which upon integration leads us to the general solution for the Eq. ([2]), for 
the choice < 4A, of the form 



x{t) = Coe-2"*cos(a;(t-(5)), (31) 
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where Co (> 0) is the second integration constant, as it should be for the 
under damped case. We also note that on substitution of the solution (!3T!) 
into ([20]), one gets I = H = (llogCo^ - ^aS). 



4.2 Case: 2 > AX (r ^0,1,2 cases): 

In this case we have the following Hamilton's equation of motion 



1 (r-1) 

X = (p) (!-'■) - ^ '-ax. (32) 

r 

(r-l) 

p = -ap 33 

r 

Again it is quite straightforward to check that Eqs. (1521) and fl55]) are equiv- 
alent to for a2 > 4A (r 7^ 0, 1, 2). 

Integrating the above equations, we obtain the general solution as 



x{t) = Coe-^°* + C^ . \. e'^^\ (34) 

(r — 2)a 

where Cq and Ci (> 0) are integration constants. Eq. (IMI) is the well known 
solution for the overdamped case. Again, we check that from (120|) that / = 

H=('--r)aCoCt''^. 



4.3 Case: 3 = AX (r = 2): 
Here the Hamilton's equations are 

la .a 
X = - - -X, P= -P- (35) 
p 2 2 

The general solution in this case is 

2:(t) = (Co + Cif)e-^°*, (36) 

where Cq and Ci (> 0) are integration constants. Here I = H — 1 = 1 — 
iaCoCfi -logCi. 
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44 Case: 4 a 0, A = (r = l): 



From (l25i) . the canonical equations for this case are 

X = e^, p = —a, (37) 
leading to the correct solution for this case, 



xm = Co - -e-"(*-*°\ (38) 
a 

where Cq and to are integration constants. Here I = H = oCq. 

Also it may be that in each of the four cases for the associated canonical 
equations, the divergence of the flow function 



A = |^ + |^=0. (39) 

OX op 

as it should be in view of the existence of conserved Hamiltonian function. 



5 Standard Hamiltonian form 



The Hamiltonian forms given in Eq. fl25|) for the various parametric ranges 
appear rather complicated: however, they can be brought to 'standard' forms 
by effecting appropriate canonical transformations. 



5.1 Case: 1 < 4A.- 



With the canonical transformation P = ^\og{x'^) and U = —xp, the Hamilto- 
nian 
into 



2 

nian for the present case, H = | log[x^sec^(u;a;p)] — can be transformed 



H = P+^ log[sec\ujU)] + ^U. (40) 
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5.2 Case: 2 (a^ > AX) and Case: 3 (o? = AX): 



For these two cases, the Hamiltonian functions given in Eq. (!25|) get trans- 
formed through the canonical transformation x = —Pe~^ and p = as 



5.3 Case: 4 a ^ 0, A = (r = 1): 

In this case, the Hamiltonian H = + ax transformed into the form 



through the canonical transformation x = ^ and p = —all. 

Since the canonical variables are now separated, the new Hamiltonians can be 
utilized for applications including the quantization. 

Before proceeding further, we wish to point out that an explicit time indepen- 
dent integral for the damped harmonic oscillator ([2]) was given by Lopez and 
Lopez [23,24]. However, no explicit Hamiltonian could be associated with it. 
We trace the reason for this in Appendix-A. 



6 Quantum description 

Having obtained appropriate Hamiltonian formalism, we now proceed to iden- 
tify a possible quantum description of the damped linear harmonic oscillator. 
For this purpose, we first rewrite the classical Hamiltonian in a symmetrized 
form so as to ensure hermiticity when the dynamical variables are represented 
by linear operators acting on a quantum Hilbert space. Thus we write 




> 4A (r ^ 0,1,2), 
= 4A. 



(41) 



H = P + e 



(42) 
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H 



^ \og[x'^ sec^ {xp + px))] — ^{xp + px), o? < 4A 

- ir-lla{xp + px), 
log(p) — ^{xp + px), = 4A (r = 2 

ef + ax, A = (r = 1) 

+ a = (r = 0) 



> 4A (r ^ 0,1,2), 
2 _ ^ \ _ o^ (43) 



In the above x and p are now considered as linear operators satisfying the 
cannonical commutation rule (CR) 



[x,p] = iH. 



(44) 



Obviously, one can look for a coordinate representation in which x is treated 
as a variable and p as the operator —ifi^ so that the CR (jH]) is satisfied 
or in a momentum representation one may treat p as a c-number and x as 
the operator x = iH-^ satisfying (jH]). For the present problem we choose a 
momentum representation and write the time-dependent Schrodinger equation 
as 



(45) 



Next we look for the stationary states 



i£t 

'i'{p,t) = tp{p)e fi 



(46) 



so that (H5|) reduces to the time-independent Schrodinger equation in the 
momentum space: 



(47) 



Here £ is identified as the eigenvalue corresponding to the Hamiltonian oper- 
ator H. In the following we first consider the pure damped case, that is, A = 
and then the remaining cases (except the case < 4A). We wish to note one 
can consider the quantization of the transformed Hamiltonians (H0l) - (l42l) also 
which is discussed in Appendix B. 
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6.1 Case:lX = 0: 



In the case A = 0, the Hamiltonian H = + ax is quantized so that the time 
independent Schrodinger equation in the momentum space becomes 



H'i/j = £^^ (ian4- + e^ip = Sip. 
dp 



(4J 



Solving the above spectral problem for any real value of£^,— oo<£^<oo, we 
obtain the wavefunction in the form 



'0(p) = C'lexp 



ah 



(49) 



where Ci is the normalization constant. From equations ( l46l) and ( 149|) we get 
the plane wave like solution 



'^{p,t) = Ci exp 
as the final solution of (1471) . 



i£t ^.feP-pS 
H \ ah 



— oo < S < oo 



(50) 



6.2 Case: 2 a^ = AX: 



For the case = 4A, the eigen value equation becomes 



Hilj = Si)^ (|log(p) - v^{l + 2p^))^ = Si) 



(51) 



Introducing the transformation of the form -0 = ^/je'fta'-^"^^^-*-'^ into the above 
equation we get 



di/j 
^ dp 



1 2£- 

2^^ ha. 



-oo < £^ < oo. 



(52) 



The above eigenvalue equation is exactly of the same form as that of the 
eigenvalue equation in the coordinate space for the toy model for quantum 
damping introduced by Chruscinski [25], see Eq. (6.2), in a rigged Hilbert 
space. 



Consequently the generalized eigenfunctions can be given as 
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i'ip) — 



1 



-tP± 



.28 1 

'aH~ 2 



where following Ref. [25] the tempered-distributions are defined as 



(53) 



p^, p> 0, 
0, p < 0, 



(54) 



P- 



0, p>0, 
\p\\ P<0, 



with A e C. The resultant wavefunction becomes 



(55) 



P± 



.2£_ 
ah 



(log(p))^ 
ha 



(56) 



In Ref. [25], Chruscinski has shown that the complex eigenvalues of the un- 
derlying Hamiltonian correspond to the poles of energy eigenvectors when 
continued to the complex eigenvalue plane and that the corresponding gener- 
alized eigenvectors may be interpreted as resonant states and these resonant 
states are responsible for the irreversible quantum dynamics. Consequently, 
similar arguments can hold good here also. Detailed description will be given 
elsewhere. 



6.3 Case: 3 a^ > 4A; 



In this case we have the following eigen value equation to be solved. 



Hip ^ Si; 



fr - 1) , , (r-2) (r - 1) ^, 



d 



(r-2) 



2r -Ml + 2p^))^ = ^^- (57) 



As we did earlier let us introduce the transformation 



■0 = tpe 



r{r — 1) 
" ah{r - 2)2^""' 



(58) 
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into equation (1571) and we obtain the eigenvalue equation 



^ dp 



rS 



11 
+ - 



[l-r)aH 2 



-oo < £^ < oo. 



(59) 



Correspondingly, again following Chruscinski [25] the generalized eigenfunc- 
tions can be given as 



2(1 — r)7raH 



P± 



rS 



1 



'l-r)ah 2 



(60) 



Then the wavefunction becomes 



^(p) 



2(1 — r)7raH 



rS 



l-r)ah 2 

p± ^ e 



1\ Y r(r-l) ^ 



ah{r - 2)' 



;P' 



■ (61) 



Again the properties of the generalized eigenfunctions discussed in Ref. [25] 
hold good here also. 



Finally for the case < 4A in Eq. (H3l) . it is not yet clear how to solve the 
linear eigenvalue problem. However, another possible approach is to make use 
of the canonically transformed Hamiltonian of the form PUI) as indicated in 
Appendix B. 



7 Conclusions 



In this paper, we have identified explicit time independent first integrals for 
the damped linear harmonic oscillator valid in different parameter regimes 
using the modified Prelle- Singer approach. We have constructed the appro- 
priate Lagrangians and explicit Hamiltonians from the time independent first 
integrals and transformed the corresponding Hamiltonian forms to standard 
Hamiltonian forms using suitable canonical transformations. We have also 
indicated that a quantized description can be developed in the momentum 
representation using the Hamiltonian forms and map onto known quantum 
mechanical toy models of damped systems. We expect that further analysis 
of the quantum case can reveal much information on the quantum damped 
oscillator. 
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A Appendix 



A set of first integrals for the damped harmonic oscillator ([2]) was reported in 
Refs. [23,24] in the form 



K = -{x^ + axx + Xx^)e-''^. (A.l) 
2 



where 



G 



ft an 



2x+ax 



log 



(a-Va^ -4\)x+2x 
{a+Va''~4\)x+2x 



Ax 

. 2ax+ct^x ' 



< 4A (underdamped) 
> AX (overdamped) (A. 2) 

^ — 4A (critically damped) 



a 



Note that in Refs. [23,24] the value of G for the critically damped case is given 
as G = but the correct form is G 



All the first integrals given in Eq. (lA.ll) are related to our first integrals 
(Eq. (l20ll -underdamped. overdamped, critically damped cases) in the following 
manner (after using suitable transformation formulas) 



-lrV4A-a^1 



(2-'-) 
2 



\ log(2K) - ., tan 
H}(2i^) 
-|log(2ir), 

where r is the quantity defined in Eq. (fT9l) . 



< 4A (underdamped) 
> 4A (overdamped) (A. 3) 
a"^ = 4A (critically damped) 



In Refs. [23,24], the authors have found the Lagrangian and an implicit Hamil- 
tonian only for underdamped case (a^ < 4A) and concluded that the dissipative 
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harmonic oscillator has nonexplicit Hamiltonian formulation. For the other 
cases it was stated that it is not possible to get a close expression for the 
Lagrangian. However in the present work we have found the Lagrangian and 
exphcit Hamiltonian for all the three cases. 

The reason for not getting the Lagrangian and generalized linear momentum 
from the first integral flA.ip reported in Refs. [23, 24] is as follows: 

For example, the first integral / given in Eq. (!20l) for the overdamped case 
is related to the first integral K reported in Refs. [23, 24] through the rela- 
tion / = {2K)^^^ . Our first integral / for the overdamped case unam- 
biguously leads to the generalized linear momentum, Lagrangian and explicit 
Hamiltonian from Eqs. fl?I]) and (1251) whereas integral K does not. The reason 
for the above is not difficult to understand. For example, let us consider the 
first integral I = x"^ + uP'x^ for the linear harmonic oscillator which gives the 
conjugate momentum p = x and Lagrangian L = x^ — u'^x'^ (vide Eqs. (pT!) 
and fl23l) ). However, / = [x"^ + w^x^)", where n is some arbitrary parameter, 
is also a first integral for SHO. By substituting this first integral / into the 
expressions (1211) and fl23|) . the linear momentum and Lagrangian cannot be 
obtained in terms of any closed form. In order to find the Lagrangian and 
Hamiltonian from the known first integral, one has to take appropriate form 
of the first integral. Our PS method gives directly the correct form of first 
integral needed to find the Lagrangian and Hamiltonian of the system. 



B Quantum description in standard Hamiltonian form 

Let us consider the Eqs. (HOj) - (l42l) . which are all of the form 

H = P + giU), (B.l) 
where 



9iU) = 



l\og[sec^{ujU)] + ^U, a2>4A 



e-°^, A = 0. 



T^e—r, >4A (r ^0,1,2), 

(r—2)a ' \ / 1 1 / 1 



The Hamiltonian (]B.1|) is now quantized in the new coordinate space, so that 
''du 



here P = —ifi-Mj and the corresponding time independent Schrodinger equa- 



tion becomes 
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h%Ij = EiIj^ {-ifi^ + giu))i) = Si). 



Introducing the transformation -0 = ifje 



J 9(U)d 



into fIB.SP we get 



fB.3) 



Solving the above one can get the form of ip, namely, 



(B.4) 



Then the wavefunction ■0 becomes 



m) = Cie 



Here 



(B.5) 



(B.6) 



g{U)dU 



^{Uhg{sec^{ujU)) + 2f/log(l + e^*^^) 

+ ff/2-^(u;f/2 + iL^2(-e2*^^))), < 4A 

> 4A 

A = 



Q(r- 2)2^ ' 



(B.7) 



where Lin{z) is polylogarithm function [26]. In order to relate this form of 
the wavefunctions in terms of the original variable one has to make use of the 
appropriate canonical transformations. 
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